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Abstract. This paper addresses the three-dimensional Navier-Stokes equations for an in- 
compressible fluid whose density is permitted to be inhomogeneous. We establish a theorem 
of global existence and uniqueness of strong solutions for initial data with small H 5-norm, 
which also satisfies a natural compatibility condition. A key point of the theorem is that 
the initial density need not be strictly positive. 
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1. Introduction 

In a number of important applications to fluid mechanics, such as in geophysical fluid 
dynamics, the Navier-Stokes equations are called upon to describe situations in which a 
fluid is inhomogeneous with respect to density, however it is essentially incompressible to a 
great extent. Such cases occur within the outer core and mantle of the earth, for example, or 
else within the atmosphere in a regime of flows for which dynamic effects of compressibility 
do not play a principle role. Ignoring any possible buoyancy or other body forces, flows of 
this character satisfy the following system of inhomogeneous incompressible Navier-Stokes 
equations; 

{d t p + div(pw) = , 
d t (pu) + div(pw ® u) - pAu + VP = , 
div u — . 

It is a standard mathematical model to pose the initial and boundary value problem, giving 
the following data; 

/ 12 ) (p,u)\ t =o = (po,u ) in ft , 

u = on <9f2 , u(x, ■) — > as \x\ — >■ oo . 

Here p, u and P denote the density, velocity and pressure of the fluid respectively, and p is 
the positive viscosity, which is assumed to be a constant. In this paper, we consider ft is 
either a smooth bounded domain in R 3 or the whole space M. 3 . 

For the initial density with positive lower bound, the inhomogeneous equations (II. II) have 
been studied in the sixties and seventies by the Russian school, see [IJ[2TJ[2I| and the many 
references that they contain. It was proved that a unique strong solution exists locally for 
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arbitrary initial data. Moreover, these papers also establish global wellposedness results 
for small solutions in dimension N > 3, while for the two dimensional case they establish 
the existence of large strong solutions. More recently, there have been many subsequent 
contributions to this theory, obtaining global wellposedness for initial data belonging to 
certain scale invariant spaces, see for example [IJ[T3j[T8]. 

On the other hand, for initial data which permits regions of vacuum, i.e. regions where 
the density p vanishes on some set, the problem becomes much more involved. Authors 
including DiPerna and Lions [T5|l25] prove the global existence of weak solutions to (11. ip in 
any space dimension, see also [T4"|I2"T] and the references in their work. As pointed in [12J, 
the major difficulty for existence of strong solutions which admit regions of vacuum is the 
lack of an appropriate estimate for d t u, since d t u only appears in the momentum equation 
with a possibly degenerate coefficient p. One way to avoid this is to estimate Vd t u instead, 
for which one pays the price of being required to impose a compatibility condition. In fact, 
Choe-Kim [12] prove the following theorem, 

Theorem 1.1. Assume that the initial data (po,«o) satisfies 

(1.3) < po G Li DH 2 , uq G -Dq' 2 fl D 2 ' 2 , 
and the compatibility condition 

(1.4) pAu - VP = \fPo9o , div m = in Q , 

for some (Po,go) G D 1,2 x L 2 . Then there exists a time T > and a unique strong solution 
(p,u,P) to the initial boundary value problem 111. ip - 111. 2ft satisfying 

p G L°°(0, T; L§ fl H 2 ), u G C([0, T] ; Dq 2 n D 2 > 2 ) n L 2 (0, T; D 3 > 2 ) , 
d t u G L 2 (0, T; Dg) , and Jpd t u G L°°(0, T; L 2 ) . 

It is also interesting to study the regularity criterion for such strong solutions, see [TT1I22] 
and its references. In particular, Kim [22] proved that if T* is the blowup time of a local 
strong solution, then necessarily 

f T * 2 3 

(1.5) / ||w(s)||r,9,°o ds = oo , for any (p,q) , - + - = l,3<g<oo, 
Jo V q 

where L q,oc is the usual weak Lebesgue space. This is effectively a criterion of regularity, and 
is similar to the Serrin criterion [26] for a priori regularity condition for the homogeneous 
Navier-Stokes equations. 

However, like the usual homogeneous Navier-Stokes equations, the question of global exis- 
tence of strong solutions remains open for large initial data. A partial result due to Kim [22] 
asserts that if H^oUd 1 - 2 is small enough, then a strong solution exists globally in time and is 
unique. In this paper, we generalize this result to the case of data within a class of critical 
Sobolev spaces. 

There is a scaling invariance of the system (II. ip . namely if (p, u, P) is a solution associated 
with the initial data (po,«o); then 

(p A , u x ,P x ) = (p(\\\x) ,A M (A 2 t,Ax),A 2 P(A 2 t,Aa;)) 
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is a solution to p.ip - (ll.2p associated with the initial data (po(Xx), Xuo(Xx)). A scale invari- 
ant function space H* is one for which || (p A , U\) \\h* = \\(Pi, u i)\\h*- For example, scaling 
invariance implies that 

IMI ■ i = \\u\\\ . i , 

II 11^ II A||^ 2 , 

which is our principle example. Our main result in this paper is that if ||wo||^i is small 
enough, then a strong solution exists globally in time. More precisely, 

Theorem 1.2. Assume the conditions in Theorem \l.l\ Then there exists a small constant 
e depending on p = \\pq\\l°°, P and the domain Q, such that if 



then the unique local strong solution constructed in Theorem li.il exists globally in time. 
Moreover, this global strong solution satisfies the following decay properties: 



The global result above relies on the smallness of the scale invariant norm of the initial 
data. In terms of the well known result of Fujita & Kato [16] for the homogeneous Navier- 
Stokes equations, our result can be considered as its generalization to the inhomogeneous 
case. The homogeneous Navier-Stokes equations have been a much more active topic of 
inquiry in this direction, see [H1E3II1E2] and their references. The key ideas leading to these 
results come from a variety of refined estimates of Stokes semigroup combined with the 
contraction principle for a Picard iteration scheme. However, in the present work one can 
not directly apply the same methods, due to the presence of the additional factor p before 
the term d t u. Instead we pursue a different strategy of proof, based on the method of energy 
estimates, and using the parabolic property of the system of equations. The main idea goes 
as follows: we first assume that the scale-invariant quantity HVuH^^ is less than an a priori 
constant bound 2; then using this assumption we then prove, under the smallness assumption 
on the initial data, that in fact the quantity is less than 1. Since ||Vu||i4 L 2 is initially 0, 
which is less than 1, then it remains less than 1. On the other hand, the boundedness of 
||V , u|| i 4 L 2 implies the local solution can be extended, according to some a priori estimates. 

The proof in this paper is inspired by a similar result for 3D compressible Navier-Stokes 
equations due to Huang, Li & Xin [19], in which they proved the global wellposedness of 
classical solutions with small initial H 2 +£ -norm. 

Remarks: Compared to the global existence results in [lj lT3~l[18] . for which the initial density 
is required to be a small perturbation of a constant, our result does not require an assumption 
of smallness on the density variations, and in fact it allows for the presence of regions of 
vacuum. 

In the case of periodic boundary conditions on a torus T d , the same result can be proved 
if the intial data has zero momentum, following the strategy as for the case of a bounded 
domain. Indeed, the same conclusion holds for classical solutions if the initial data are regular 
enough and satisfy some higher-order compatibility conditions. The proof of this statement 
is almost the same except that some additional high-order estimates are required. These 






Vw|| L 2 < Mr* , 
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estimates have been given in [22] . In the case in which f2 is a smooth bounded domain in 
IR 2 , it is already known that global strong solutions exist without a smallness assumption; 
this was proved first for the case without vacuum [I] and subsequently for the case with 
vacuum [20] . 

Our result holds in particular in the case of a constant density p = 1. It therefore gives 
a new proof of the well known Fujita-Kato theorem [16J , which uses only energy estimate 
methods in the argument. 

Theorem 1.3. Assume uq 6 if 1 (R 3 ) and div uq = 0. There exists a positive constants, such 
that if || Wo || a-^ < £ > then there is a global strong solution to the homogeneous Navier-Stokes 
equations with the property that 

M GL7([0,oo);^ 1 )nLL(0,oo;^ 2 ) . 

This article is organized as follows: section 2 sets the notation, and contains several 
definitions and basic lemmas. In section 3 we give the proof for the case Q = R 3 , while the 
proof for the case of a bounded domain is presented in section 4. 



2. Notation and interpolation lemmas 

The homogeneous and inhomogeneous Sobolev spaces are defined in the standard way. 
For 1 < r < oo and fcGN, 

U = L r (tt), 

D k > r = {u e L} oc : \\V k u\\ L r<oo}, \\u\\ D k,r = \\V k u\\ L r , 

W k,r = L r p D k,r ? jjk = yykj ? jjKr = ^ c \ osme in t he norm of D k ' r . 

12 i 

In particular when f2 is a bounded domain then D ' = Hq. Set 

C% a = {u e C °° : div u = 0} , 

and set D^J 1 = C^ a , with the closure taken in the norm of D 1,2 . The fractional-order 

homogeneous Sobolev space H S (R 3 ) are defined as the space of tempered distributions u 
over M. 3 for which the Fourier transform Tu belongs to L i 1 oc (IR 3 ) and which satisfy 

H^ (R 3 ) :=y" 3 ierwoi 2 de<oo. 

And H S {£L) is the restriction of H S (K 3 ) to the domain f2. We refer to H S (Q) by H s , inde- 
pendent of whether Q is the whole space or not. 

Some well known interpolation results which relate Lorentz spaces and the classical Sobolev 
spaces are presented, which can be found in |8,28,29j. For every < 9 < 1, 1 < q < oo, and 
normed spaces X ,Xi, we use the notation [Xo,Xi]e >q to denote the real interpolation space 
between X and X\. 

Lemma 2.1. Let s, s , Sj el. If s = (1 — 9)s + 9s\ with some < 9 < 1, then 



[H S0 ,H S ^ 2 = H S . 
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And the following inclusion relation holds, if 1 < qi < q2 < oo, then 

[H s \H s % qi Q[H s \H s % m . 
Lemma 2.2. Let < 9 < 1, 1 < po,Px,p < oo, 1 < q , qi, q < oo, and ^ = + then 

(2.8) [L Po ' 9o ,L Pl ' 91 ] 0i? = L p ' q . 

This identification of the interpolation spaces (\2.8\i is also valid if L Pl,Ql is replaced by by 
L°° (L Po,qo , L pi ' qi and L p,q are Lorentz spaces, whose definition can be found in JEj). When 
q = p, L p,p is the Sobolev L p space. 

The above interpolation lemmas clearly generalize to vector-valued Lorentz spaces. For 
example, 

Lemma 2.3. Let < 9 < 1, 1 < po,Pi < oo, 1 < qo, q%, q < oo, and - = + Then 

(2.9) [L po ' 90 (0,T;L 2 ),L Pl ' 9l (0,T;L 2 )] ei9 = L M (0,T;L 2 ) . 

The statement (EH]) is also valid after replacing L Pl ' qi (0, T; L 2 ) by L°°(0,T;L 2 ). 

Bounds for linear operators on interpolation spaces are the concern of the next lemma. 

Lemma 2.4. Suppose that E , Ei, F Q , Fi are Banach spaces and O<0<1, l<g< oo ; 
and that A is a linear operator from Eq + E\ into Fq + F\. If A maps Eq into Fq, with 
— -^oll^lUo f or Q tt w e -^0; and maps E\ into F\, with \\Aw\\f 1 < LiHiwH^ for all 
w £ Ei. Then A is a continous linear operator from [Eo,Ex]g iq into [Fo,Fi\g tq and one has 

PHI[Fo,F!] e , 9 < Ll" 6 L{\\w\\ [E ^ El]eq . 

In this paper, the constants C, C%, C2, C3, C4, C5, • ■ ■ may depend upon /i,p, Q and as 
conventional in analysis, they may change from line to line. 

3. Proof of Theorem 11.21 for = M 3 

In this section, we focus on the case Q = M. 3 . 

3.1. A Priori Estimates. This subsection establishes several a priori estimates for strong 
solutions to the Cauchy problem, which will play a key role in extending local strong solutions 
to global ones. 

Given a strong solution (p, u, P) on M 3 x [0, T], define 

(3.10) A(t) = \\Vu\\ LH0>t . L2) , < t < T . 

Theorem 3.1. Under the assumptions of Theorem there exists a positive constant e 
depending only on jj and p, such that if \\uq\\^i < e, and (p,u,P) is a strong solution to 
(HHJ -dOIl, satisfying 

(3.11) A(T) < 2, 
then it in fact holds that 

A(T) < 1. 

The remainder of this subsection consists in proving this key result. 
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Lemma 3.2. Let (p,u,P) be a strong solution of fll.ip - fll.2p . Then for all <t <T, 
(3.12) = IIpo||l°° = P ■ 

Proof. The mass equation is in fact a transport equation, owing to the fact that div(-u) = 0, 
from which (I3.12p follows. □ 

Lemma 3.3. Let (p,u,P) be a strong solution of flLip -f TO]) . Suppose that 

A(T) < 2, 

then there exists some constant C\ depending on p, p, such that 
(3-13) ||V«||i4 (0)T . L 2) < CilKH^i , 

and it holds that 

(3.14) sup £||V«(t)||| 2 <C p \u \ 2 dx. 

te[o,T] J 

Note that C\ does not depend on T . 

Proof. First, let's consider the following linear Cauchy problem for (w,P), 

pd t w — pAw + (pu ■ V)w + VP = , 

(3.15) { div w = , 

w(x, 0) = Wq(x) . 



Suppose that Wq satisfies the conditions assumed for uq as prescribed in Theorem ll.il then the 
existence and uniqueness of strong solution to (13. 15[) has been proved in |12j . Straightforward 
energy estimates tell that 

(3.16) ~ J p\w(T)\ 2 dx + p J \\Vw\\ 2 L idt<~ J p \w \ 2 dx <C{p)\\w \\ 2 L 2 . 
Multiplying (I3.15p i by d t w and integrating over IR 3 , one gets by Sobolev embedding that 

J p\d t w\ 2 dx + J \Vw\ 2 dx = — J " (pu ■ V)w • w t dx 

(3.17) <C(p)\\u\\ L e\\Vw\\ L s\\^d t w\\ L2 

<^ J p\d t w\ 2 dx + C\\ VM||^2||Vw|| L 2||Vw|| i 6 . 

Notice that the momentum equation can be written as 

— pAw + VP = —pd t w — (pu ■ V)w, 

where the left handside is viewed as the Helmholtz-Weyl decomposition of the right one. 
From this equation, 

p\\Aw\\ L 2 < \\pd t w\\ L 2 + \\(pu ■ V)w\\ L 2. 
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It follows from Calderon-Zygmund inequality and the Sobolev embedding theorem that 

||Vw|| L 6 < l7||V 2 w|| L 2 

(3.18) <C\\pd t w\\ L 2 + C\\{pu ■ V)w\\ L 2 

<C\\pd t w\\ L 2 + C\\Vu\\ L 2\\Vw\\l4Vw\\\ 6l 
which, together with Young's inequality, implies that 

(3.19) l|Vw|| L6 < C\\V 2 w\\ L 2 < C\\pd t w\\ L 2 + C\\ Vw|| 2 2 || Vw\\&. 
Inserting (13.191) into (13.171) . we get 

(3.20) l - J p\d t w\ 2 dx + ^ J \Vw\ 2 dx<C\\Vu\\\ 2 \\Vw\\\ 2 . 
By Gronwall's inequality and the assumption A{T) < 2, 

(3.21) J j p\d t w(s)\ 2 dxds + || Vtu(t) Uia < Ce c \\Vw \\ 2 L2 . 

For fixed (p,u), the map from wq to Vw(i) is linear. Furthermore, by Lemmas 12.11 and 
E3i: 4 (0,T;L 2 ) = [L 2 (0,T;L 2 ), L°°(0, T; L 2 )} i 4 and C U 2 ,^ 1 ! , then one gets, upon 

2' 4 

combining the estimate (13.16!) and (13.211) . and using Lemma 12.41 that 

( 3 -22) ||Vw|| L 4( 0) T;i2) < Cl||Wo||^l . 

Notice that G\ does not depend on T, since one can scale [0, T] to [0, 1]. Consequently, 

||Vm||l4( 0jT;L 2) < Ci||« ||^i. 
Regarding the decay of ||Vrx||x,2, upon multiplying (I3.20p by t, one gets that 

t J p\d t u\ 2 dx + — J t\Wu\ 2 dx< J \Vu\ 2 dx + Ct\\Vu\\tf\\Vuf L 2. 

If A(T) < 2, then by Gronwall's inequality, 



(3.23) 



T 

t\\y/pdtU\\] ? dt+ SUp t\\ Vu(t) 

o te[o,T] 

T 



<C7^ \\Vu\\ 2 L 2 dt < C J p \u \ 2 dx, 
where the basic energy inequality for u is utilized. 



□ 



Proof of Theorem \3.1[ The conclusion of Theorem 13.11 will follow if we let e < 1/Ci, which 
is given in Lemma [3.31 □ 

Remark 3.1. According to the proof of Lemma \3. 3\ the H 2 -norm of Uq can be replaced by 
. 1 

B2 A -norm, which will make the main result more refined. 
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From this point on, we will consider only the small data problem, assuming that the initial 
data satisfies the condition ||w || 1 < £, as in Theorem 13. II The notation C is used to denote 
a positive constant, which may depend on T and the initial data, and it may change from 
line to line. The computation is standard, with respect to that in [22], but we sketch it here 
for completeness. From the proof of Lemma 13.31 we know that 



(3.24) 



/ Ha/p^Hl 2 dt + mII^ u IIl~(o,T;I, 2 ) ^ C|| Vw ||i2 < C . 
Jo 



Lemma 3.4 (Estimates for ||*/]o^u||£,2 and || Vw||#i). Under the assumption of Theorem \3.1[ 

we have 

(3.25) sup [\\^pd t u\\ 2 L2 + \\Vu\\ 2 m ] + f \\Vd t u\\ 2 L2 dt < C . 

te[o,T] Jo 

Proof. Differentiating the momentum equation with respect to t, multiplying by dtu, and 
then integrating over M 3 , one can obtain that 



p\d t u\ 2 dx + p J \Vd t u\ 2 dx 



1 d 
2di 

(3.26) 

= — 2 / pu - Vd t u ■ d t u dx — I pu ■ V(u • Vw ■ d t u) dx — I (pd t u ■ V)w • d t u dx 



It follows from Sobolev embedding theorem and Gagliardo-Nirenberg inequality that 

— 2 / pu - \7d t u ■ d t u dx 
(3 ' 27) <C||Vw|U 2 ||V^||| 2 ||Vp^||| 2 

<^llv^ni 2 + L7||v M || 4 L2 ||yp^lli 2 - 

By Sobolev embedding and the estimate (13. 19[) , the second term can be estimated, 
— J pu ■ V(u ■ Vu • d t u) dx 

</H«l|v M n ail< |^+//>l»l 2 |v 2 »l|aH^+/*l 2 |v«llvaHrfx 
^3,28) <c|klUe|l VM IU 2 ll Vw IU 6 ll^ M IU 6 + c|I m II16||v 2 m||l2|I^ m IU 6 

+ C7||u||i«||Vu|Ue||V^u|| za 
<L7||V M || 4 L2 ||V 2 M ||i 2 + ^||V^||i 2 

<C\\Vu\\U\ P d t u\\h + L7||Vn||i 2 + |||V^||| 2 . 
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For the third term on the right handside of (13. 26p . utilizing Gagliardo-Nirenberg inequality 
and (|33SD , 



/< 



(pd t u ■ V)u ■ d t udx < C\\yfpd t u\\ L 2\\Vu\\ L 6\\y/pd t u\\ L 3 
( 3 - 29 ) <C\\^pd t u\\\ 2 \\Vd t u\\% (\\pd t u\\ L 2 + \\Vu\\ 3 L2 ) 

10 n 

<c\\^pd t u\\i, + c\\jpd t u\\l4Vu\\ A L * + ^iiv^ni, . 

o 

Collecting all the estimates (I3.27l) - (l3.29p . one gets that 



(3.30) 



lit I P^ u \ 2dx + ^ I \VdM z dx 



it 
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< C\\Vu\\U\y/pd tU \\l 2 + C\\^pd t u\\l 2 \\yfpd t u\\ 2 L2 + C||V 

Utilizing Gronwall's inequality and the estimate (13.241) . we obtain that 

(3.31) sup \\^pd t u\\ 2 L2 + [ T \\Vd t u\\ 2 L2 dt < C. 

te[o,T] Jo 

According to the estimate ( 13.191) . ( 13.31[) implies that 

(3.32) sup ||V 2 m|| L 2<l7 sup || yfpd t u\\ L 2 + C||Vu||| 2 < C . 

te[o,T] te[o,T] 

It completes the proof. □ 
Lemma 3.5 (Estimate for ||Vp||#i). Under the assumptions of Theorem \3.1\ we know that 

(3.33) sup \\Vp\\m + / II V 3 n||| 2 dt < C . 

*6[0,T] Jo 

Proof. Differentiating the mass equation with respect to Xj, j = 1,2,3, we find that djp 
satisfies 

(3.34) d t djp + u ■ Vdjp = -dju ■ Vp. 

Then multiplying (13.341) by djp, integrating over M 3 , and summing over the index j, one gets 
that 

(3.35) ^_ J \Wp\ 2 dx <C J \Vu\ ■ \Wp\ 2 dx < C\\ Vu\\ L oo\\ Vp||| 2 . 

To derive the appropriate bound for J Q T || Vii||L°° dt, we make use of the elliptic estimates 
related with the momentum equation. In fact, 

II V 2 m|| L 4 < C\\pd t u\\ L , + C\\(pu ■ VH| L 4 

<C\\Vd t u\\l4^d t u\\l 2 + C\\u\\ L i2\\Vu\\ L s 
< C\\Vd t u\\l4^d t u\\l 2 + C\\Vu\\l4Vu\\l 6 , 
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which together with Lemma 13.41 gives that 

(3.36) ||V 2 M || L2 (0,T;L4) < C . 

By virtue of the Gagliardo-Nirenberg inequality, 

||Vm||l« < l7||Vm||| 6 ||V 2 m||| 4 , 

hence HVwHii^T;!, 00 ) < C. Consequently, by Gronwall's inequality, (I3.35[) gives a bound for 
1 1 Vp 1 1 ip.. Similar argument shows that 



(3.37) 



j t J \V 2 p\ 2 dx<C J (|Vw||V 2 p| 2 + |V 2 w||Vp||V 2 p|) dx 



< C||VM|| L ~||V 2 p|li2 +L7||V 2 M || L 3||Vp|| i6 ||V 2 p||L 2 

< l7(||Vm|| L oo + ||V 2 m|| l3 ) ||V 2 p||i 2 . 



Note that ||V 2 m|| L 3 < ||V 2 «||^ 3 ||V 2 «||J1 3 . Combining this fact with (13321) and flQj) , we get 
that 

V 2 u e L^CTjL 3 ) . 

Hence, the Gronwall inequality gives a bound for ||V 2 p||l2. 

Finally, using the reguality theory for Stokes equations, one can obtain that 

II V 3 u|| L2 < C{\\pd t u\\ H i + \\pu ■ Vu\\ m ) 

(3.38) 

< C (|| Vp|| L 3 + 1) (|| Vd t u\\ L * + || Vu||ffi) , 
which implies that || V 3 m||l2( 0jT;L 2) < C. □ 

3.2. Proof of Theoerem 11.21 With the a priori estimates in subsection 3.1 in hand, we 
are prepared for the proof of Thorem 11.21 

Proof. According to Theorem II. 1| there exists a T* > such that the Cauchy problem 
(jl.ip - (|1.2[) has a unique local strong solution (p,u,P) on R 3 x (0, T*], where T* depends on 
llpoll^l H 2i ||Vtto||ifi and H^Hl 2 - We will show that this local solution extends to a global 
one. 

It follows from the integrability property of the local strong solution that A(0) = 0. Hence 
there exists a T\ G (0, T*) such that (13.1 ip holds for T = T\. Set 

T* = sup{T| (p, u, P) is a strong solution on M 3 x (0, T] and A(T) < 2} . 

Then T* > T x > 0. 

The claim is that T* = oo, for otherwise T* < oo, which we will assume for an argument 
by contradiction. First, it follows from Lemmas 13.41 and 13.51 that 

||V«(t)||Hi < C7(T*) , ||p(t)|| L 3 nif2 < C(T*) , for every < t < T* , 

where C(T*) depends on T* and the initial data. Secondly, 

pAu — VP = pd t u — (pu ■ V)u, 



INHOMOGENEOUS INCOMPRESSIBLE NAVIER-STOKES EQUATIONS 



11 



which implies that 

<\\sfpd t u\\ L 2 + \\(y/pu ■ V)u\\ L 2 

^llV^wll^ + C^llVtillli-H^lli 
<C{T*) . 

These uniform estimates allow us to construct a new unique local strong solution from 
(p(T 2 ), u(T 2 )) for any < T 2 < T*, for which the life span is uniformly bounded from below 
by some AT. By uniqueness, on the interval [T 2 , T*), this new strong solution coincides with 
the original one. Choose some T 2 close enough to T*, such that T 2 + AT > T*, giving a local 
strong solution (p, u, P) on (0, T 2 + AT]. 

Next, we show that (13. lip holds on (0,T 2 + AT]. For the strong solution (p,u,P) on 
[0,T 2 + AT], let 

(3.39) T** = sup {T| A(T) < 2} . 
It follows from Theorem 13.11 that in fact 

A(T**) < 1 . 

Hence T** = T 2 + AT and (13. lip holds on (0,T 2 + AT]. In fact this is a contradiction to 
the assumption that T* < oo is the maximal time, since T 2 + AT > T*. Furthermore the 
decay property of the solution is implied in the proof of Lemma 13.31 Hence, the proof is 
complete. □ 

4. Proof of Theorem 11.21 for the case of a bounded domain 

In this section, we assume that the domain Q is bounded and smooth domain in M 3 . As 
before, for strong solutions (p, w,P) to (Il.ip -( T1~2|) on Q x [0,T], define 

(4.40) A{t) = ||V«||£4 ( o )t . ^ , 0<t<T. 

Proceeding as in the above section, we establish a priori estimates for (p, u) which guar- 
antee the extension of the local solution. 

Theorem 4.1. Under the assumptions of Theorem \l.l\ there exists some positive constant 
e depending only on p, Q and p, such that if ||wo||^.i < s and (p,u,P) is a strong solution 
to flUU-flOP, satisfying 

(4.41) A(T) < 2, 
then it in fact holds that 

A{T) < 1. 

The proof for Theorem 14.11 is similar to that for Theorem 13. 1[ with some slight changes due 
to the presence of several additional lower order terms that appear because of the existence 
of a boundary. These terms are controlled using the exponential decay of the kinetic energy. 

Lemma 4.2. Let (p,u,P) be a strong solution to (ll.ip - (ll.2l) . Then for every < t < T, 

IIpWIU- = IIpoIU- = P ■ 



p~2 (pAu - VP) 
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Lemma 4.3. [Modified Energy Estimate] Let (p,u,P) be a strong solution to (jl.ip - fll.2D . 
Then it satisfies 

(4.42) \\Vpu(t)\\v < \\^u \\ L 2ex V {-C(fi,n,p)t} . 
Proof. According to the basic energy estimate, 

(4.43) ~ J p\u\ 2 dx + p J \Vu\ 2 dx = . 
By the Poincare inequality, 

f P \u\ 2 dx < c{n,p)\\Vuf L 2 , 



which, when incorporated into (14.431) . yields 

d f 12 2u j i i . . 

dtJ pH dx + cWJ) lM dx -°- 

Consequently 

p(t)\u(t)\ 2 dx < I po\uo\ 2 dx ■ exp{—C(p,il, p)t} . 



Lemma 4.4. Let (p,u,P) be a strong solution to (jl.ip ~ (jl.2p . Then 

(4.44) sup t / p\u\ 2 dx+ t |Vm| 2 dxdt < C / p \u \ 2 dx . 

te[o,T] J Jo J J 

Proof. Multiplying the momentum equation by tu and integrating over Q, 

t d 



2dt 

which implies that 



p\u\ dx + pt I \Vu\ dx = , 



sup - p\u\ 2 dx + p / t \Vu\ 2 dxdt = - j j p\u\ 2 dxdt 

te[o,T] 2 J Jo J 2 , 

Taking (14.421) into account, one finds that 

(4.45) sup t / p\u\ 2 dx + p / t | Vu\ 2 dxdt < C I pq\uq\ 2 dx , 



□ 



te[o,T] 

which completes the proof. □ 
Lemma 4.5. Let (p,u,P) be a strong solution of the equations (ll.ip - (ll.2l) . Suppose that 

A(T) < 2, 

then there exists some constant C\ depending on p, p, Q, such that 

( 4 - 46 ) ||Vtt||L4(0,T;L2) < Ci||mq||^J ) 
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sup £||V-u(i)||f 2 < C I pol^ol 2 ^- 
efo,Tl ' J 



and it holds that 
(4.47) 

te[o,T] 

Note that C\ does not depend on T . 

Proof. The proof is similar as in the case of the whole space. We consider the following linear 

system for (w,P), 



(4.48) 



pd t w - pAw + (fw ■ V)w + VP = , in Q x (0, T] 
div w = , infix [0, T] , 
w = , on dQ x [0, T] , 
w(x, 0) = wq(x) , in fl . 



Suppose wq satisfies the conditions assumed of uq in Theorem 11.11 then the global existence 
and uniqueness of strong solution to (I4.48[) is known. Straightforward energy estimates tell 
that 



(4.49) 



1 



p\w(T)\ 2 dx + /i / \\Vw\\ 2 L2 dt < - / Pq\wq\ 2 dx <C(p)\\wo\\ 



|L2 



Multiplying (I4.48p i by d t w and integrating over Q, one gets by the Sobolev embedding 
theorem that 



p\d t w\ 2 dx + 



p d 



Vw| 2 dx — — / pu ■ Vw ■ d t w dx 



(4.50) 



2 dt 

<C\\u\\ Le \\Vw\\ L s\\^d t w\\ L 2 



<- A J p\d t w\ 2 dx + C\\ Vu||| 2 || Vit?||i2 1| V 



W\\ H i . 



From the regularity theory for stationary Stokes system, it follows that 
||Vw|| H i < C(\\pd t w\\ L 2 + \\(pu ■ V)w\\ L 2) + C\\Vw\\ L 2 
< C\\pd t w\\ L 2 + (7||Vu|| Z 2||Vw||t 2 ||VHIIri + C\\Vw\\ L 2 

which implies 

(4.51) \\Vw\\ m <C\\pd t w\\ L 2 + C\\Vu\\ 2 L 4Vw\\ L 2 + C\\Vw\\ L 2 . 

Inserting the estimate (14.511) in (I4.50p . then 

J p\d t w\ 2 dx + /i— J \Vw\ 2 dx 

< C||Vu||ia||Vlw|||a + C1|Vw||i 2 ||Vw|| 2 2 

< C?||Vu||ta||Vlw|||a + HIVHIia • 



(4.52) 
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Utilizing Gronwall's inequality and Poincare inequality. 



(4.53) 



f \\^d t w\\l 2 dt + f i\\Vw(T)\\l 2 <Ce c \\Vw \\l 2 +C l 2 [ \\Vw\\ 2 L2 dt 
o Jo 

<Ce c \\Vw \\ 2 L2 +C J Po\w \ 2 dx<Ce c \\Vw \\ 2 L2 . 



By interpolation, L 4 (0,T;L 2 ) = [L 2 (0, T; L 2 ), L°°(0, T; L 2 )} 1>4 and #2 C [L 2 ,^ 1 ]^, so 
that 

( 4 - 54 ) I|Vw||l4 (0iT . £2) < CiHiWoll^ , 

which also holds for u. 

Furthermore, multiplying (I4.52p by t we have 



t J p\d t u\ 2 dx + /i-^ J t\Vu\ 2 dx<p J \Vu\ 2 dx + Ct||Vtt||| 2 + £tt||Vu 
By Gronwall's inequality and Lemma I4.4[ for every t G [0,T], 

J s j p\d t u(s)\ 2 dxds + t\\Vu\\ 2 L2 < C J (1 + s)\\Vu\\ 2 L2 ds <C j 
which completes the proof. 



|2 
\L 2 



p \u \ 2 dx, 



□ 



The proof of Theorem 14.11 is thus complete, if we set e < \jC\. Finally, the proof of 
Theorem 11.21 for bounded domain case follows the same logic as the analog result in section 
3, we omit the details for reasons of brevity. 
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